Equilibrium under flow 
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Using information theory we derive a thermodynamics for systems evolving under a collective 
motion, i.e. under a time-odd constraint. An illustration within the Lattice gas Model is given 
for two model cases: a collision between two complex particles leading to a incomplete relaxation 
of the incoming momentum, and a self-similar expansion. A semi-quantitative connection with the 
determination of thermodynamical quantities in multifragmentation reactions is done showing that 
they are affected in a sizeable way only when the flow dominates the global energetic. 
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Heavy ion collisions represent a unique opportunity to 
probe the interdisciplinary field of phase transitions in 
mesoscopic systems for which a non standard thermody- 
namics giving rise to negative heat capacities is predicted 
W . An indication of negative heat capacities has been ex- 
perimentally found in multifragmentation reactions. 
These results can be considered as a first step towards a 
quantitative determination of the nuclear phase diagram. 
Negative heat capacities have also been observed in sim- 
ulations of self-gravitating systems Q and in the melting 
and boiling of clusters Q . 

An important conceptual problem linked with multi- 
fragmentation experiments is that the outcomes of a nu- 
clear collision are not confined in an external container 
but dynamically deexcite in the vacuum. In the absence 
of boundary conditions, collective flows can be present at 
the fragmentation time. Because of this time odd com- 
ponent, one may doubt about the applicability of equilib- 
rium concepts. However, we will show that information 
theory allows a thermodynamically consistent description 
of open finite systems in evolution under a collective flow. 
We will apply this formalism to the Lattice Gas Model || 
for the particular cases of a memory of the entrance chan- 
nel (transparency) and radial expansion. In both cases 
the properties of the system appear to be affected in a 
sizeable way when flow dominates the global energetics. 

In the Gibbs formulation of statistical mechanics an 
equilibrium is defined as a collection of different mi- 
crostates all corresponding to the same macrostate (or 
statistical ensemble) defined through the average value 
of (a collection of) collective variables. A practical real- 
ization of a Gibbs ensemble is given by the collection of 
different snapshots of an isolated ergodic system evolv- 
ing in time, provided that the observation time is much 
longer than a typical equilibration time . This type of 
statistical ensemble is not very useful for short-lived open 
systems as the transient excited states formed in a nu- 
clear collisions. However, ergodicity is not the unique 
way to produce a Gibbs ensemble . Indeed, within infor- 
mation theory, an equilibrium corresponds to any ensem- 
ble of states that maximizes the entropy in a given space 
under the constraint of (a number of) observables known 
in average^). Therefore, if the nuclear dynamics is suf- 
ficiently sensitive to the initial conditions, the ensemble 



of outcomes of similarly prepared nuclear collisions can 
be considered as a statistical ensemble for which the im- 
portant observables are controlled by the dynamics (and 
by the event sorting performed on specific observables). 
As an example, in the standard freeze out hypothesis 
the configurations are fixed when the nuclear interaction 
among prefragments becomes negligible; then the ensem- 
ble of events is characterized by several variables such as 
the energy and the spatial extension R 2 which for open 
systems is fluctuating event by event. This can be ac- 
counted for by considering the ensemble average < R 2 > 
as a state variable and introducing a Lagrange multiplier 
A closely related to a pressure. In this information the- 
ory approach, the state variables are determined by the 
dynamics and can also be time odd quantities such as 
transparency or radial flow. 

Let us first consider a symmetric head-on collision with 
a too short reaction time to fully relax the incoming mo- 
mentum. This situation seems to be verified at relativis- 
tic energies |Q. It corresponds to the observation of an 
additional one body state variable, the memory of the ini- 
tial momenta < ep z > , where p z is the momentum along 
the beam axis and e = — 1(+1) for the particles initially 
belonging to the target (projectile). Let us assume that 
the total energy E is also known only in average. The 
maximization of the entropy leads to the partition sum of 
the incomplete momentum relaxation ensemble (IMRE) 
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where the indice i stands for the i-th particle while (n) 
counts the events. (3 and a are Lagrange multipliers as- 
sociated to the constraint of the degree of incomplete 
stopping < ep z > and of the total center of mass energy 
< E > . The relative probability of an event (n) results 
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where Uij is the two body interaction and we have in- 
troduced po = ma/(3u z . The average kinetic energy is 
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< E k > I A = 3/(2/3) +pl/(2m) while the equation of 
states related to a leads to < ep z >= Ap . In the limit 
Po = 0, these expressions correspond to the usual canon- 
ical ensemble while in the general case they can be inter- 
preted as two thermalized sources with a non zero relative 
velocity 2po / m along the beam axis u z . In actual heavy 
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FIG. 1: Fragment size distributions in the IMR Lattice Gas 
model. Upper part: constant temperature plus an increasing 
proportion of directed linear momentum. Lower part: micro- 
canonical calculations. Full lines: canonical and microcanon- 
ical calculations with an isotropic momentum distribution. 

ion experiments the centrality selection criteria imply a 
sorting of data according to the total deposited energy ||] 
or to variables which are strongly correlated to it. This 
means that a total energy conservation has to be imple- 
mented to eq.(||). In this case the equations of state are 
not analytical but can still be numerically evaluated. 

This theory can be extended to take into account other 
observables. if for example the average momentum dis- 
persion along the beam axis is known, an extra term 
dPiz^ 2 can be added in eq.(Q) to control the collective flow 
fluctuation, while a term XR^ 2 can be used to impose an 
average freeze out volume. For simplicity in the follow- 
ing calculations the longitudinal momentum dispersion 
has been kept fixed, (f}/2m + <5) _1 = 0.04e. 

To understand the effect of transparency on the eval- 
uation of thermodynamical quantities we have used the 
Lattice Gas hamiltonian Q where occupied sites on a 
three dimensional cubic lattice interact via a constant 
closest neighbors coupling e. This simple but numeri- 
cally solvable model is isomorphous to the Ising model 
in the grancanonical ensemble and constitutes therefore 
a paradigm of standard equilibrium statistical mechanics 
with first and second order phase transitions. Calcula- 
tions in the IMRE eq.(|l]) are performed for a system of 
A = 216 particles at a subcritical pressure A = 3.3 • 10~ 4 



and clusters are defined within the standard Coniglio- 
Klein prescription 0. Since in this classical approach 
the configurational and kinetic partition sums are fac- 
torized, for a given [3 the lattice configurations will be 
independent of the degree of transparency pq. However 
the active bond probability will explicitly depend on pa 
meaning that cluster observables can be affected by the 
incomplete relaxation. Figure la shows the cluster size 
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FIG. 2: Average partial energies (left) and variances (right) 
in the microcanonical ensemble (full circles) and in the IMRE 
(open circles) at fixed energy. 

distributions for a temperature = 0.65e which cor- 
responds to the transition temperature in the canonical 
ensemble, and different quadrupolar deformations in mo- 
mentum space 

Ap = V2< \ Pz \ > - < \p±\ > 
P 2(V2< [p,| > + < \p x \ >) 

A quantitative comparison with experimental data would 
require to fix the Lagrange parameters A, /?, po from each 
specific set of data and to include quantum effects and 
symmetry as well as Coulomb terms in the hamiltonian, 
however from this illustrative example we can clearly see 
that partitions are affected by a collective longitudinal 
component and a higher degree of fragmentation does 
not necessarily imply higher temperatures but can also 
be consistent with an increased degree of transparency 
of the collision. To quantify this statement, Figure lb 
compares the clusters size distributions of the IMRE at a 
fixed total energy with the standard microcanonical en- 
semble (fixed energy, spherical momentum distribution) 
at the same energy. It is clear that thermal agitation 
is much more effective than transparency to break up 
the system: at the energy where the microcanonical en- 
semble predicts a complete vaporization of the system, 
a residue persists if the non relaxed momentum compo- 
nent is as large as the relaxed one. This is in qualita- 
tive agreement with the trend observed in central colli- 
sions at relativistic energies On the other side up 
to 10% transparency the distributions are almost identi- 
cal, meaning that when the velocity difference between 
the quasi-projectile and the quasi-target is of this order, 
the debate on equilibrium based on the number of emis- 



3 



sion sources jl4] is an academic question. The important 
influence of the collective component on the cluster dis- 
tributions means that, if the incoming momentum is not 
completely relaxed, the methods used to determine ther- 
modynamical quantities can be strongly biased. Indeed 
if energy is divided into a kinetic et and an interaction e% 
component, the average kinetic energy can be used as a 
microcanonical thermometer while the normalized partial 
energy fluctuations are linked to the microcanonical heat 
capacity Q . The interaction energy can be estimated as 
the Q- value of the detected fragments . Figure 2 shows 
the deformation induced by transparency on the first and 
second order moments of partial energies that are used 
to infer the temperature and heat capacity in the micro- 
canonical sorting of fragmentation data j| . The Q- value 
can be computed in the Lattice Gas model in a liquid 
drop approximation as 
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where M^> is the multiplicity of the n-th event and a^, 
a s are the volume and surface energy coefficients. The 
average kinetic energy and the average Q-value obtained 
with eq. @) are shown in the left part of figure 2 for the 
microcanonical ensemble and for a system with the same 
total deposited energy and an increasing degree of trans- 
parency Ap/p. The energy range explored corresponds to 
the phase transition region. The increased probability of 
heavy clusters with an increasing degree of transparency 
shown in figure 1 induces an overestimation of the average 
kinetic energy, hence of the estimated temperature. The 
peak in the partial energy fluctuation which signals the 
phase transition j^] disappears in this extreme scenario 
of an increase of energy going entirely into the relative 
motion of the two sources, and the normalized fluctua- 
tions are systematically lower than in the microcanonical 
ensemble where the incoming momentum is fully relaxed. 
Since a negative heat capacity corresponds to abnormally 
high partial energy fluctuations ||] , Figure 2 implies that 
an incomplete relaxation of the incoming momentum can 
prevent the observation of negative heat capacity. 

Another important form of collective motion in heavy 
ion collisions is radial flow which starts to be observed in 
central collisions around 30 A.MeV incident energy) 14 
and becomes the dominant fraction of the detected en- 
ergy in the relativistic domain E|. We can describe 
this dynamical situation in the framework of informa- 
tion theory as an equilibrium with non random directions 
for velocities which are preferentially oriented in the ra- 
dial direction. In the canonical formulation this corre- 
sponds to two independent observations of the average 
energy < E > and a the average local radial momentum 
< p r (r) > . The probability of a microstate (n) reads 
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FIG. 3: Upper left: effective pressure (full line) and aver- 
age volume (dashed line) as a function of the collective radial 
velocity. Upper right and lower part: fragment size distri- 
butions in the expanding Lattice Gas model. Distributions 
without flow (full black lines) are compared with distributions 
with 10% (dashed black), 30% (full grey) and 60% (dashed 
grey) contribution of radial flow at the same temperature (up- 
per right) and at the same total energy (low part). Symbols: 
calculations with flow at the same thermal energy and average 
volume as the full black lines. 



where r, is the center of mass position of the i-th par- 
ticle and (3, 7 (r) are Lagrange multipliers. Imposing in 
the local equation of state < p r >— dlogZ/dj that the 
observed velocity is self-similar < p r (r) >= mar we ob- 
tain j(r) = —(3ar which gives for the argument of the 
exponential in the probability (0) 
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with the local radial momentum po(r) — maru r . The 
situation is equivalent to a standard Gibbs equilibrium 
in the local expanding frame. This scenario is often 
invoked in the literature []lof to justify the treatment 
of flow as a collective radial velocity superimposed on 
thermal motion; however eq. (^contains also an addi- 
tional term oc r 2 which corresponds to an outgoing pres- 
sure. This term cannot be avoided in the selfsimilar sce- 
nario; indeed the total average energy under flow reads 
< E >= 3/(2/3)+ < E pot > +ma 2 /2 < R 2 > where we 
have used the fact that the scalar product between the 
flow and the thermal component is in average zero. The 
probability under flow eq.(|^) diverges at infinity reflect- 
ing the trivial dynamical fact that asymptotically parti- 
cles flow away. This divergence should be cured by in- 
troducing an external confining pressure which is not a 
mathematical artifact but has to be interpreted as dis- 
cussed above as a Lagrange multiplier imposing a finite 
freeze out volume. Eq.(^) has then to be augmented by 
a term — XJ2i r 1 with A > mo? jlT leading to a positive 
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pressure coefficient A e // = A — m{3a 2 /2. 

The effective pressure X e ff as we U as the associated 
average volume (normalized to the ground state volume 
Vq = A) are shown in the upper left part of figure 3 
as a function of the collective radial velocity for a given 
pressure A = 1.23 • 10~ 2 and temperature /3 _1 = 0.65 e. 
The Lagrange parameter A e // being a decreasing func- 
tion of a, the flow reduces the effective pressure so that 
the critical point is moved towards higher pressures in 
the presence of flow §. However one can see that the 
effect is very small up to ma ~ .6 (which corresponds 
to as 40% contribution of flow to the kinetic energy). In 
this regime the cluster size distributions displayed in the 
upper right part of figure 3 are only slightly affected. On 
the other side if the collective flow overcomes a thresh- 
old value AEk/Ek ~ 50% the average volume (dashed 
line in fig. 3a) shows an exponential increase and the out- 
going flow pressure leads to a complete fragmentation of 
the system (dashed grey line in fig. 3b). Again an oriented 
motion is systematically less effective than a random one 
to break up the system. This is shown in the lower part 
of figure 3 which compares for a given A distributions 
with and without radial flow at the same total deposited 
energy: for any value of radial flow equilibrium in the 
standard microcanonical ensemble corresponds to more 
fragmented configurations. 

Concerning heavy ion collisions, it is important to 
stress that the pressure A as well as the other state vari- 
ables are consequences of the dynamics. They cannot 
be accessed by a statistical treatment but have to be 
extracted from simulations and/or directly infered from 
data itself Different models assume that fragmenta- 
tion occurs in an average freeze out volume which may 
depend on the thermal energy but does not depend on 
flow. This is true if the system fragments at the turning 
point of its expansion (A e // = 0) p5| or when the in- 
teraction between fragment surfaces becomes negligible 
or more generally insufficient to modify the N-body 
correlations O]. In this case the presence of flow does 
not affect the configuration space and can only modify 



the partitions because of the modified microscopic bonds 
among particles taken into account by the Coniglio-Klein 
algorithm. However this effect is negligible as already ob- 
served by Das Gupta et al.JHI and shown in the lower 
part of Figure 3. In this figure the symbols represent the 
size distributions with collective flow at the same thermal 
energy and average volume as the standard microcanon- 
ical results (full lines). Even for the largest amount of 
flow considered the two distributions are identical, mean- 
ing that in this hypothesis all thermodynamical analysis 
of fragmentation data stay valid in the presence of even 
strong collective flows . 

In this paper we have shown that information theory 
allows to take into account generic collective motions 
in a rigorous statistical way leading to the formulation 
of a consistent "out of (static) equilibrium" thermody- 
namics. We have particularized this general statement 
to two kinds of collective motion which are particularly 
relevant in nuclear multifragmentation, namely partial 
transparency and radial selfsimilar flow. We have ap- 
plied the formalism to the Lattice Gas model which has 
already shown some pertinence Jl6| to the multifragmen- 
tation phenomenology. We have shown that the disor- 
dered thermal motion is always more efficient than the 
collective motion to break up the system. A consequence 
of that is that normalized partial energy fluctuations for 
a system under flow are always reduced respect to the 
standard microcanonical expectation^ and the negative 
heat capacity signal experimentally observed || cannot 
be due to the possible presence of a collective component 
neglected in the data analysis. Transparency leaves the 
configurational energy unmodified, while radial flow, be- 
cause of the coupling between r and p space, is shown 
to be equivalent to an external pressure leading to an 
increased value of the critical point 0. Concerning the 
extraction of thermodynamical variables from fragmen- 
tation data, radial flow can be simply subtracted from 
the detected energy |H| while transparency can be ne- 
glected only if the relaxation of the incoming momentum 
exceeds about 90%. 
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